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We extend the fermionic functional renormalization group (fRG) flow for symmetry-broken states 
to a case of competing ordering tendencies in a one-dimensional model with attractive interactions. 
First, as an introduction to the fRG method with symmetry-breaking counterterms we study the 
reduced BCS model. Then the ground state phase diagram of an one-dimensional model with 
competing orders is analyzed by fRG and compared with bosonization. We present a self-consistency 
argument which allows us to find the optimal counterterm for general models that cannot be solved 
exactly by mean-field theory. We show that the fRG is capable of reproducing the quantum phase 
transition between dominant pairing correlations and CDW order in a meaningful way. 



I. INTRODUCTION 



Many of the phenomena in condensed matter physics can be understood as consequences of electronic correlations. 
Superconductivity, magnetism, and the quantum Hall effects are just a few examples of how collective behavior of 
electrons leads to macroscopic effects. In many cases the main effects can be ascribed to some kind of longer-range 
order. It is thus an important task of condensed matter theory to shed light on the question why and how electrons 
order in different ways. If the order parameter is known or restricted to a few possibilities, much progress can been 
made within mean- field approaches like the BCS theory of superconductivity. On the other hand, the question if and 
how the electrons order in microscopic models such as the Hubbard model is still not resolved in general. One reason 
for this question still being open is that various ordering tendencies strongly influence and compete with each other. 

The onset of long-range order is usually connected with a phase transition. One of the most powerful tools for the 
investigation of phase transitions is the renormalization group (RG), laid out in its most common form by Wilson in 
197tf4 Wilson's RG is formulated in terms of the partition function of the system, that is, a functional integral over 
the degrees of freedom of the system. The idea is to integrate out the degrees of freedom one after another depending 
on their wavelength, in contrast to perturbation theory, where all modes are treated simultaneously. One begins at 
high energies and short wavelengths and works one's way down to low energies, taking into account the modes that 
are integrated out through changes in the parameters of the action of the remaining modes. This process is called 
renormalization, and the change of the parameters from high to low energies is called the RG flow. 

Later, Polchinski 2 formulated the renormalization group as an exact flow equation for the partition function of 
the theory. In the following, this was generalized to other generating functionals, such as the one for one-particle 
irreducible (1PI) vertices that will be used in this work. The flow equations for this generating functional were 
first derived for scalar fields by WettericfP, and later for fermionic systems^HE ! By n0 w it has been widely used 
to investigate many-fermion systems in various contexts (e.g. Refs. I7|8|9|10|ll|121 In these systems the dominant 
question is which type of symmetry-breaking takes place at the transition, while the precise nature of the transition is 
often only of secondary interest. For the RG schemes used in this context, the term functional renormalization group 
applies to the method in its present form in a twofold way: one uses a flow equation for a generating functional of 
fermionic fields to obtain flow equations for vertices, which are in turn functions of the generalized quantum numbers 
(wavevectors, Wannier orbitals) of the degrees of freedom. This has to be contrasted with usual RGs which are in 
general flow equations for just a small number of coupling constants. 

The functional RG was in particular employed to study instabilities of the Fermi surface of the 2D Hubbard model 
towards different symmetry broken states such as antiferromagnetism, superconductivity, and ferromagnetism by 
various researchers (e.g. Refs. I7|8|9|1Q|1T]) . However, for technical reasons, it was not possible in these studies to 
obtain quantitative results for physical quantities like the order parameter, as the flows used were restricted to the 
symmetric phase. Instabilities in these studies would only show up in the form of diverging susceptibilities. In 2004, 
Salmhofer et al. 13 found a way to continue the RG flow into phases with broken symmetry by including a small 
external field in the model, which is explained in detail in section \\IB As these external fields bias the results, they 



are by no means an optimal method for calculating order parameters. However, Gersch et alP^ extended this idea by 
replacing the unphysical external field by a counterterm, an external field that is added to the interacting Hamiltonian 
but at the same time subtracted from the free Hamiltonian, so that - at least in principle - the calculation is not 



biased any more (as long as calculations are performed exactly). This will be discussed in section III 



Until now, the method of counterterms has only been applied to a model where the flow equations can be solved 
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exactl}ff^. The goal here is to see how the various approximations involved when dealing with non-trivial systems 
affect the method. More concretely, we want to investigate if the counterterm scheme a) is able to deal with situations 
where mean-field theory is not exact and b) in situations in which instabilities occur in more that one channel at 
a time. We consider the attractive Tomonaga-Luttinger model at half-filling as a simple toy model. This model 
allows numerically cheap computations, whilst at the same time exhibiting interesting behaviour with two competing 
instabilities in the region of the phase diagram considered. A further advantage is that very much is known about 
the Tomonaga-Luttinger model (as will be discussed in section IV), so that the data obtained by the functional RG 
approach can be compared with existing information. 

Of course, as the Tomonaga-Luttinger model is already well-understood, no new physical results will be presented. 
Our interest here is purely methodological, and is thought to serve as a preliminary, but important test of the method 
before applying it to more interesting situations, such as the repulsive 2D Hubbard model. Another test in the 
attractive 2D Hubbard model produced promising results that this approach can be used to get beyond standard 
mean-field calculations^. 

Before embarking onto the derivation of the necessary fermionic functional RG equations, let us comment on an 
alternative promising approach for flows into symmetry-broken regimes in interacting many-fermion systems. Instead 
of following the four-fermions interactions down to low scales, it is also possible to decouple the fermions via an adaptive 
Hubbard- Stratonovitch transformation in the relevant bosonic channels (which channels one should actually choose 
can be read off from the fermionic flows in the symmetric phase}^. This has the advantage that the large fermionic 
interactions are replaced by nearly massless bosons and can cause no further problems in the truncated flow equations. 
Furthermore collective fluctuations of the order parameter (i.e. the ordering boson) can be treated more easily. Often 
the connection to standard models of statistical mechanics such as 0(A^)-models is feasible. On the back side, the 
generation of couplings in channels not present in the initial interaction (such as d-wave pairing in the 2D Hubbard 
model on the square lattice) and the removal of regenerated fermionic interactions are technically challenging 1 ^, and 
the bosonic propagators and the Yukawa couplings to the fermions need to be approximated strongly for practical 
calculations. This in turn reintroduces some degree of bias and complicates the application to more realistic models. 
While these problems can certainly be taken care of to some degree, we feel that in addition one should try to develop a 
complementary, purely fermionic approach that is able to capture the most relevant corrections beyond the mean-field 
studies and that could be used to determine the basic gap structures of unconventional superconductors in a relatively 
straightforward setup. 



II. THE FUNCTIONAL RENORMALIZATION GROUP FOR THE 1PI VERTICES 



A. Derivation of the Flow Equations 

In this section the fRG flow equations for the one-particle irreducible (1PI) vertex functions are described. The 
derivation presented here follows Ref. 4. For alternative derivations, see, e.g., Ref. 6. The fRG is based on the 
evaluation of a generating functional describing a quantum many-body system by introducing a parameter x into 
the functional such that for some value Xi the functional can be evaluated exactly and for some other value Xf 
the functional describes the system one is interested in. This parameter is called the flow parameter. There are 
many different ways of impleneting this idea, which correspond to different choices of generating functionals and flow 
parameters. Evaluating the functional is then done by virtue of the identity 

lim f(x) = f(Xi)+ [ X f(x)dx. (1) 

X^Xf J x% 

The limit in the above formula will turn out to be important, as generating functionals in the thermodynamic limit 
may possess singularities, so the result in general depends on how Xf is approached. This corresponds to the fact 
that in the thermodynamic limit phase transitions may occur. 

The starting point for the derivation is the partition function - the generating functional of correlation functions. 
For fermionic systems it is given by the Grassmann functional integral 

where we have have also inserted a normalization factor and Grassmannian source fields. The action is now separated 
into a free (quadratic) part and an interacting part, 



sm = \u>,Qii>)+v[ii>], 



(3) 
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where the i/j are Grassmannian fields, which are viewed as Grassmann algebra- valued vectors indexed by the multiindex 
X = (x,r, c, s, . . .). The charge index c is introduced to write both -0 and ijj as components of one field. Q is the 
inverse free propagator of the system, which is a matrix in the vector space generated by X. With the introduction 
of the charge index for the fields, Q is antisymmetric. The notation (• , • ) denotes the standard inner product of the 
vector space generated by the fields ipx- The normalization factor Zq with these conventions is given by the integral 
over the free action only, Zq = y/detQ. The square root comes from using the charge index. 

Next we introduce the real flow parameter \ into the free action, Q — > Q x , where the exact form of the parameter 
dependence of Q x need not be specified for the derivation. This induces a flow parameter dependence of the partition 
function as well, Z — > Z x . In order to arrive at the desired flow equations for the generating functional of 1PI vertices, 
which is also called the effective action, we have to perform a series of transformations on Z x . 

First, we define the generating functional of connected correlation functions by 




We define the field expectation values as 



dW x r i 

-ftfM- (5) 



This naturally induces the inverse r] [</>] = r] x [4>\. Equipped with these definitions one can perform the Legendre 
transformation of W x [rj] with respect to 77 to arrive at the effective action 

r x [<t>] = -w x m}-(<P,vl<P})- (6) 

In order to derive a flow equation for T x [</>], we differentiate eq.Qwith respect to x to find 

t x = -d x W-(f,,d v W)-(<p,f,), (7) 

where all dependences on <\> and 77 have been omitted and the dot denotes the total derivative with respect to %• O n 
inserting the definition of (j) given in EqjEJ one finds that all terms except the first one in Eq. [7] cancel, so that 

t x [<P] = -d x W x [r l [4>]\. (8) 

Hence, it is clear that it is sufficient to evaluate the derivative of W x , which is readily done by taking into account 
the definition of W x in eq.Q. Writing Z x [rf\ = exp{W x [r]]} on has 

d x Z x = d x e w * = (d x W x )Z x (9) 



and (omitting some subscripts x) 

d x z = Vd^Q(d x ^=)z- 1 - /Jp 



e -sw+(i),(/i) 



= {^( a *^)-5(^*)} z - 

The first term still has to be evaluated, which is readily done employing the well known identity from linear algebra, 
In det = Tr In, to obtain 

d Y . = d Y exp I — -Tr In Q 1 

= -l^vm- (10) 
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Finally, one has to employ the definition of Z = e w into eq.(10) and evaluate the derivatives acting on Z in the second 
term of eq.(10), which leads to 



Z = 



((« QS V W) + Tr (QS^W)) Z, 



(11) 



where S^S^W is understood as a matrix in the vector space indexed by the usual multiindex, where the indices of the 
sources r\ define the matrix indices. Putting all together, we have arrived at 



d x W = - 1 - {Tr (QQ- 1 ) + (s n W, Q5 v W^j + Tr (Q5 v 5 v w} } 



(12) 



which now has to be translated into a formula for f by replacing (rj) by ((/>, <p) and W by T. This is done using the 
identities (see, e.g. Ref. [18]) 



fyr = 77, 



Upon using them, one finds that 



f = ^Tr (QQ- 1 ) + \ (4>, Q4>) + ^Tr (Qf" 1 ) , 



(13) 



(14) 



where f^- denotes the matrix d^.d^T. This is the central result of this section, obtained first in Ref. 4 . An equivalent 
formula for the scalar case was worked out before in Ref. |3l The trace in the last term can be understood as a one- 
loop diagram, i.e. the flow can be evaluated quite simply for cases where one has a good ansatz for T. For fermions 
however, T needs to be expressed in monomials of Grassmann fields which leads to a infinite set of coupled one-loop 
equations for the vertex functions. This will be explained in the remainder of this section. 
The expansion of T into a basis of monomials in the fields <j) is done by writing 



m>0 



with 
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(m)r 



ml 



d m X 7m (X)0 m (X). 



(15) 



(16) 



Here X = (Xl, • • • ,X m ) is an m-tuple of multiindices, and m (X) = <fi(Xi) • • • <j>{X m ). The 7 m (X) are chosen to be 
totally antisymmetric, reflecting the fermionic nature of the fields involved. Furthermore, we similarly expand 



8<j>{X) 5<f>(Y) 



m>2 



where 



j (m \X,Y;4>) = ^^yy j d m X' 7m (*, Y, X')<F" 2 (X') 



(m 

One finds in particular that 7^) does not depend on </>, so that 

s 2 m 



5<P(X)S<P(Y) 



72 



(17) 



(18) 



(19) 



4>(X)=<t>(Y)=0 



From eq.(13) it follows thus that 



72 = Q~\ 



(20) 
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where Q denotes the full propagator, so that by defining 

s 2 m 



54>{X)5<p{Y) 



= 12 {X,Y) + T{X,Y;, 



72 



(i + W)) . 



With this definitions, the flow equation, eq.(14) reads 



(21) 



(22) 



This equation is still non-polynomial because of the last term. This term, however, may be written as a power series 
as 



Tr 



QQ (l + Qt (0)) 1 = Tr (GO) - Tr (QQGT (0)) 

^(-i) p Tr [gq (gm) 1 



p>2 

By defining the so-called single scale propagator S as 

S = -GQG, 



(23) 



(24) 



and comparing equal powers of the fields in eq.(22) using eq.(23) we finally arrive at the flow equations for the 1PI 
vertices: 



7 (6) W>) = 2 Tr 



Tr 
Tr 



S^\4>)Gi {i) (4>) 



S {^Gl {i) + 7 (4) <?7 (6) 



-Tr 



5^(4)^(4)^(4) 



(25) 
(26) 
(27) 



Actually, one has an infinite hierarchy of ordinary differential equations, but as in actual computations only the lowest 
equations are used, the equations for the higher order vertices are not displayed here. 

In order to compare coefficients the RHS of the flow equations have to be antisymmetrized, because the LHS are. 
Carrying out this antisymmetrization one arrives at the following set of equations for 72 and 74: 

(28) 
(29) 

(30) 
(31) 



i 2 {X u X 2 ) = Q{x 1 ,X 2 ) + \ J &X 3 AX i S{X i ,Xz) li {X 1 ,X 2 ,X 3 ,X A ), 
7 4 (X) = \j dY 1 dY 2l6 (X,Y 1 ,Y 2 )S(Y 2 ,Y 1 )-^ J d 4 YL(YB(X, Y), 



where Y = (Yi, • • • , Y4) and similarly for X, and furthermore 



l(y) = 5(r 1 ,r 2 )G(y 3 ,n) + 5(y 3 ,n)G(y 1 ,y 2 ) 

B(X,Y) = 74(X 1 ,X 2 ,r 2 ,r 3 )74(n,^i,X3,X4) 

74 (Xi , x 3 , y 2 , r 3 )74 (n, n , x 2 , x 4 ) 
+7 4 (x 1 ,x 4 ,r 2 ,r 3 )74(n,yi,x 2 ,x 3 ). 



It is thus obvious that one possibility to truncate the infinite hierarchy of flow equations is to neglect the contribution 
of 76 in eq.(29), and to solve only the closed system of equations for 7 2 and 74. The higher (non-lPI) correlation 
functions are then given by tree diagrams consisting of full propagators and the 1PI vertices given by the 74. It is 
this strategy that will be pursued in the following. A justification for dropping 7 6 under certain conditions is given 
inRef. HI 
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B. Flows into Phases with Broken Symmetry 



The main focus of this work is the calculation of order parameters corresponding to states with spontaneous 
symmetry breaking (SSB). In the setting described above, the fRG can not be used to calculate these order parameters. 
The technical reason for this is that if the initial system is invariant under some symmetry group G, then so will 
be all the 1PI vertices in the expansion of T. Now the RHS of the flow equations contain only sums of products of 
the 1PI vertices. But if all the vertices are invariant under G, then so are products of them. Thus no symmetry 
breaking terms can be generated during the flow. If the flow starts in the symmetric phase, it will remain there. In 
the case of a second order phase transition this will lead to divergences, because in this case susceptibilities diverge, 
and this divergence is reflected in the four-point correlation functions, as a susceptibility is renormalized by the four- 
point correlation function of the fermions. This method has already been used extensively to study the 2D Hubbard 
modePPPl In this case, of course, the flow cannot be performed until the end, so that the low-energy behaviour of 
the system can not be extracted precisely. 

One workaround for this problem has been proposed in Ref. [T3land consists of the application of a small external 
field that corresponds to the symmetry breaking one expects. For example, if one expects an instability towards 
5-wave superconductivity, one would add to the initial Hamiltonian an (infinitesimally small) s-wave pairing field 

A ext> 



A extE(4 T c U) +hx - 



In this approach, the symmetry is broken from the very beginning of the flow, so that 1PI vertices do exist that break 
G, thus allowing for a flow of other G-breaking terms, such as the order parameter. On the other hand, of course, 
results obtained with this method will depend on the magnitude of the applied field, and only in the limit of vanishing 
external field the desired unperturbed results are recovered. Of course this idea corresponds to the above mentioned 
fact that generating functionals in the thermodynamic limit possess singularities that mark phase transitions. The 
procedure of introducing a small external field is similar to choosing the correct ordering of the limits N — > oo and 
A — » in statistical physics. 

The model considered in Ref. [13] is the so-called reduced BCS (rBCS) model, which is described by a many-fermion 
system where the only interaction is via the Cooper channel with zero total momentum. Because there exists only 
one channel, mean-field theory is exact for this model. One can show that this implies that the truncation of the 



flow equations discussed in section [HA] is also exact. This makes this model a good testing ground for studying the 
effect of an external field without having to care about large couplings, effects of the truncation and the like. In the 
remainder we briefly summarize the results of Ref. [13] 
Concretely, the model is described by the Hamiltonian 

H = w ° - V E /® /®4i c U4r<£*P ( 32 ) 

k,p 

where f(k) describes the angular momentum of the pairing and will for simplicity be set to f(k) = 1 in the following, 
v is the interaction strength and V is the volume of the system. This choice of f(k) corresponds to s-wave pairing 
and does not lead to a loss of generality, as in the rBCS model the different pairing channels decouple. 
Next, the external field is introduced as above by adding the term 

Wext = ^E 4i c U + St E C K C -%V ( 33 ) 

k k 

describing an external s- wave pairing field E, to the Hamiltonian. The symbol E instead of the more common A is 
chosen here because the external field is treated as an initial condition for the self-energy as described below. 

In order to handle the new vertices (and lines) appearing due to the symmetry breaking, it is convenient to introduce 
the usual superconducting Nambu spinors (in the functional integral language, see Ref. [19]) und to group wavevectors 
k and Matsubara frequencies uo = ko to a multiindex k , 

*fc = ( ^ H Y = (&T . 1>-h) > (34) 



so that the matrix representation of the quadratic part of the action becomes 

Q = 



T I . (35) 



-St 
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Here Q(°) denotes the propagator in the symmetric phase. The external field is introduced as initial condition for 
the self-energy, so that E in the definition of the inverse free propagator is actually zero, and appears only in the full 
propagator, connecting fields with opposite Matsubara frequencies. 

Diagrammatically, the off-diagonal elements of the propagator can be viewed as new types of lines that cre- 
ate/annihilate two particles forming a singlet. During the flow, these new lines also lead to new vertices that do 
not conserve particle number. In the rBCS model, there are two of them corresponding to creation/annihilation of 
four particles forming two singlets. However, the anomalous terms creating particles are related to those annihilating 
particles by complex conjugation. Thus, choosing the external field to be real, their numerical values are actually 
identical. The anomalous vertex will be denoted by w in the following. As this vertex is only generated during the 
flow, its initial condition is w = 0. 

Before stating the results it is important to mention that the flow equations derived in section |II A| have to be 
modified in order to obtain the correct result for the gap. The Katanin modification^ of the truncated flow equations, 
eqns. (28) and (29), amounts to the substitution 



G 



(36) 



in the flow equation of 7^, where the dot indicates the derivative with respect to the flow parameter. The advantage 
of this modification is that it improves the fulfillment of Ward identities, which are violated in the truncated fRG 
scheme. Katanin has shown that as a consequence the flow of symmetry breaking order parameters is also improved, 
so that in the ladder (or bubble) approximation, where only particle-particle (particle-hole) diagrams are kept, the 
mean-field equation for the order parameter is recovered, which is not the case when using the original truncated 
flow equations. The Katanin modification is thus a necessary ingredient when performing flows in symmetry broken 
phases and will be always used in the remainder of this work. 

The (truncated) flow equations may be obtained from eqns. ([28]) and (29) and read 



v = B(v 2 
w = A (y 2 + w 
E = C(v + w) 

where A, B and C refer to the loop integrals given by 

A = J dk 



w z ) + 2Avw 
2A l -2Bvw 



\Q\ 2 + (Z x (s))i 

(u 2 + e 2 ) X (s) 2 



B = / dk 

c = T ^I dk \Qr- + ^x(s)r 



(|QP + (E X ( S ))2) 2 
2 X ( S ) X ( S )S|Q| 2 



(37) 
(38) 
(39) 



(40) 
(41) 
(42) 



Here we have only written the equation for the vertices with zero total incoming frequency. Vertices with nonzero 
total frequencies also flow, but more weakly and do not give rise to any self-energy corrections in the thermodynamic 
limit due to the restricted wavevector-structure of the bare interaction. We could have equally well started with 
the doubly reduced interaction where only zero total frequency pairs interact. The anomalous selfenergy remains 
frequency independent. 

Integrating these flow equations one arrives at a value of the final gap that, obviously, depends on the value of the 
external field, as shown in Fig. |IIB| on the left panel. In order to obtain good results it is thus desirable to use as 
small fields as possible. Unfortunately, there is a trade-off, shown on the right panel of Fig. II B As one lowers the 
value of the external field, the couplings grow larger. In this case, the growth of the couplings represents the massless 
Goldstone mode associated to phase fluctuations of the order parameter^. 

As the flow equations for this model are exact, large couplings (i.e. four point vertices) do not really matter. For 
more realistic problems, however, it is crucial to keep the couplings at moderate values because the truncation of the 
flow equations corresponds to a weak coupling approximation: Terms of order V 3 , where V are the couplings, are 
neglected. Thus the use of external fields puts one into the dilemma of either getting errors due to the unphysical 
field or due to large couplings. In order to meliorate the situation we will try a slightly different approach with 
counterterms. 
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FIG. 1: External field dependence of (a) order parameter and (b) couplings for the rBCS modeP^. 



III. INTERACTION FLOW AND COUNTERTERMS 



In this section we discuss the introduction of counterterms into the £RG scheme as suggested by Gersch et alP^. In 
order to do this it is necessary to abandon sharp momentum cutoffs (as will be explained below) and to use a soft or 
even flat cutoff instead. We will always use the so-called interaction flow^, so this will be topic of the next section 
before actually introducing the counterterm method. 



A. Interaction Flow 



The interaction (or flat-cutoff) flow scheme is implemented by introducing the following flow parameter dependence 
of the quadratic part of the action, Q, 

Q -> Q x = % (43) 
X 

where x is a ^-independent real number G [0, 1]. x — is the starting point of the flow. As the Green's function of 
the virtual fluctuations is Q^ 1 , this cutoff choice increases the weight of the fluctuations continuously from to 1 for 
all fc, i.e. corresponds to a completely flat cutoff. Recalling the definition of the full action, 

S X [9] = (%Q X 9) + V[9], (44) 

and assuming that V[^] is quartic in the fields, one finds by rescaling 

* = VX* (45) 
that this choice of flow parameter amounts to switching on the interactions gradually, 

S x [9] = (*,Q*) +X 2 Vm. (46) 

Thus the interaction flow scheme meets the criteria for a sensible flow scheme as (a) the partition function at x — 
can be evaluated exactly and (b) x = 1 corresponds to the physical system of interest. The reason for calling this 
scheme a cutoff scheme with an infinitely soft cutoff is that x is a constant function in momentum space and may 
hence be viewed as a cutoff function in momentum space in the limit where the slope of the cutoff is sent to zercP^. 



B. Counterterms 



The idea of introducing counterterms is simply to add zero in the form 

= a — a 



(47) 
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to the Hamiltonian (or, equivalent ly, to the action) of a quantum system. This is to obtain a different separation of 
the Hamiltonian into a 'free' and an 'interacting' part, explicitly, 

H = H +H I (48) 
Tio — > Tio + STi 
Tij — Tii — $Ti. 

For this operation to to be useful, of course the counterterm 5Ti should be quadratic in the fields. This change is very 
similar to choosing a different expansion point for perturbation theorjl^. 

Following Gersch et al. 14 , we introduce counterterms in order to break a symmetry explicitly during the flow as 
in the method using small external fields described above. The difference to the introduction of an external field is 
that with the counterterm the symmetry of the full Hamiltonian is preserved, so that no error is introduced into the 
calculation of the corresponding order parameter, hence the counterterm does not need to be small. This last point, 
however, only holds if the generating functional of the problem is evaluated exactly, because if and only if this is the 



case the two copies of STi introduced in Eq. 49 cancel. 



To see how this works in practice, the rBCS model already introduced in section II B will serve as a simple example. 
It has the additional advantage that the truncated flow equations for this model are exact, so that the result will not 
depend on the choice of counterterm (more precisely, its numerical value). In fact, the flow equations for the cutoff 



flow with an external field, eq.(51) can be taken over almost literally. To recall them, these equations are given by 

v = B (v 2 + w 2 ) + 2Avw (49) 
w A (v 2 + w 2 ) + 2Bvw (50) 
E = C(v + w) (51) 

where A : B and C refer to the loop integrals given by 

b = rV / dk i" 2 + 4Ms) 2 ( ) 



c - rV f d , Ms)x(s)nQ\ 2 m 



Here s is the flow parameter, whereas x( s ) 1S a cutoff function in momentum space. Now we use the s-dependence of 
x(s) corresponding to the interaction flow 

X(s) = 8, se[0,l]. (55) 

Next, in order to cancel the external field (or initial self-energy) E(x = 0) from the action, the counterterm A is 
added to the quadratic part Q, i.e. Q — > Q + A. Remembering that 

g- 1 = Q - E, (56) 

it is clear that this term A must be given by A = E(x = 0). Thus, in order to get the flow equations for the symmetric 
situation with counterterms from the equations with symmetry-breaking externel fields, one only has to augment the 
free inverse propagator, Q, by an anomalous term A, and the self-energy in the full propagators, which was just 
before, has to be shifted according to 

X E - X £ - A, (57) 



The expression on the right hand side is called the effective gap. These changes do not change eqns.([51|, but eqns.([54j) 
now read 

--E/^<^^) 2 

b = r V / dk (" a + ^W 8 (59) 

(|QI 2 + (£y(s)-A)2) 2 



c 



(|Q|2 + (E X (s) - A)2f 

A 3 -2 X S(A 2 + W 2 + e l) + A(( X E) 2 



^J dk |Q| 2 + (E X ( S )-A) 2 • (60) 
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FIG. 2: Flow of the effective gap for the rBCS model. The gap flows to the same value for all choices of counterterms. 



These flow equations may now be integrated, which in fact leads to a result for the gap that is independent of the 
choice of counterterm, and which is not affected by the application of an unphysical external field. In Fig. IIIB we 
show some flows for different values of the counterterm. 

We note that an alternative way to express this recipe is the following: First one adds the counterterm A to the 
quadratic part of the self-energy. Then one sets the initial condition E(% = 0) equal to A, such that in the full inverse 
propagator for x — 1 the effective gap — A is exactly cancelled. However in the flow before reaching x — 1? the 
systems appears to be in the symmetry-broken state. This allows to uncover the spontaneous symmetry breaking. 

For the reduced BCS model, the counterterm methods works just as well as the external field scheme of Ref. [T3l 
In many cases^, the counterterm method has the apparent advantage that large couplings can be avoided for a large 
part of the flow, or more precisely, by keeping the flow in the symmetry-broken state for most of the flow range, 
criticality can be avoided. This works best for systems without Goldstone bosons (e.g. discrete symmetry breakings). 
But even for breaking of continuous symmetries where the Goldstone modes require a diverging component of the 
vertex at the end of the flow, avoiding the criticality of the amplitude modes that are usually gapped at the end of the 
flow could represent an improvement. Of course, the valid question is whether also non-trivial models can be treated 
correctly with this method. In the following we will apply the counterterm flow to a model in one spatial dimension 
with competing ordering tendencies. 



IV. ATTRACTIVE ONEDIMENSIONAL MODEL - BOSONIZATION TREATMENT 



A. Model 



In general, the action for free spin- 1/2 fermions in one dimension at T = reads 



S F = ^2 J ^ x J dr^ s (^,r)(S r + e k 



■ /i)^ s O,r), 



(61) 



where e k is understood as a differential operator defined by the identification k — > id x . For low energies, E <C /i, the 
dispersion may be linearized around the Fermi points, yielding 



e k F +q ~ £k F + dk^k 



k=kp 



q = \i + v F q, 



(62) 



where vp is the Fermi velocity defined by above relation. Note also that vp at the left Fermi point differs in sign 
from vf at the right Fermi point. In order to write the effective low-energy action it is convenient to define chiral 
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components of i/j as 

R s (x,T) = Y,e lkx Mr) (63) 

fc>0 

L s (x,r) = ^e^(r) 

k<0 

where ^(r) denotes the Fourier transform of ifj(x^r) with respect to x. These new fields describe right and left 
moving electrons, respectively, and in terms of them, the action can be written as 

Sf = ^ / dxdrRl(x,r)(d T + ivFd x )R 3 (x,r) + Ll(x,r)(d T — ivFd x )L 3 (x,r) (64) 

where now vp is taken to be positive and inversion symmetry is assumed. Crossterms between R s (x,t) and L s (x,t) 
are neglected. 

Now we include interactions. As we are looking for a situation with competing orders, the attractive Tomonaga- 
Luttinger model at half- filling will be in the focus of interest. Half band filling implies that at arbitrarily small 
energies umklapp scattering is an allowed process and will be consequently included in the calculations. Traditionally, 
the interactions for this situation are classified according to the g-ology '-scheme into various channels, see Fig. [3j In 
the following, the processes are neglected, because they are not relevant for symmetry breaking. Furthermore, the 
interaction vertices in Fig. [3]do not carry spin-labels for the in- and outgoing legs. For spin-1/2 fermions, each vertex 
can in principle occur with parallel or perpendicular spins for the ingoing (outgoing) legs. However, because of the 
Pauli principle the process with parallel spins is only allowed if other quantum numbers for the fermions involved 
differ. This means that for the gs process only the perpendicular process is allowed as both ingoing (outgoing) legs are 
located at the same Fermi point. Also g^\ = — g 2 \\ holds in general so that these two vertices actually correspond to 
the same vertex. This leads to the set of vertices {#i, g 2l 9\\ , #3}, where the (i = 1,2,3) refer to the perpendicular 
process in the corresponding channel. If the couplings are invariant under spin rotations, g\\ = g 2 —g\- The interaction 
Hamiltonian with this notation (and without g±) is given by 

51 LtRl a L_ a R a - 52 LtRl a L a R_ a 

a a 

-5|| Yl LlRtL a R a + f£ (LlLl.R-.Rr + RlRl a L_ a L a ) . 




B. Bosonization 



In the bosonization procedur e 1 23 * 24 *, fermionic operators for left- and right-movers are expressed by exponentials of 
a bosonic field </>. Omitting the spin index for a moment, one writes 



-R, R 
L,L 



1 



/27ra 
1 

/27ra 



exp 



exp 



±W47T (/>(£, r) 
=fiV47T(j)(x~T) 



(65) 



For systems with spin-1/2 one needs two bosons (j) ai a G {|, I}, instead of one. However, it turns out to be more 
convenient to define instead 

(j) a = T - (66) 

P = 0t + ^. (67) 



Using the bosonization formulae, Eq. [651 the bosonized form of the full interacting Hamiltonian is readily found to 



be 



7i = 7Yo + Tii — Tt a + T~L P . 



(68) 
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FIG. 3: Classification of effective low energy interaction vertices in the Tomonaga-Luttinger model. Solid lines are located in 
the vicinity of one Fermi point (e.g. right), dashed lines in the vicinity of the other (e.g. leftj^. 



Here H p (H a ) depends on p (0 a ) only, which means that these two fields actually decouple. The expressions for the 
H pi n = p : a are given by 



dx 



dx 



Kp n 2 p + ^(d x 4> P ) 2 

J^CJ 



2tt 2 



cos 



2«* - ( V ' ST °' ; ) ' 



COS 



and the constants K p and u p are 



/2tt- (gjj T92) 
2tt + (g\\ T52) 




(69) 
(70) 



(71) 



where the upper (lower) sign refers to p = a(p = p). In particular we note for later use that for g\ < 0, K G < 1, while 
K p can be tuned to be larger than 1 by making g\\ + g2 negative. 

The Hp have the same structure, they are both of the well-known sine-Gordon form, which implies in particular that 
they exhibit a phase-transition of the Kosterlitz-Thouless type as the parameters (i.e. the couplings) are varied^. This 
phase transition can be analysed using the Renormalization Group, leading to the phase diagram for the Tomonaga- 
Luttinger model. This will be the content of the next subsection. 



C. Phase Diagram 



The derivation of the one- loop RG equations for the sine- Gordon models of Eq. 68 can be found for example in 
Ref. [24j As spin and charge sector decouple, there are two independent sets of flow equations, which in the charge 
sector are given by 

The equations for the spin sector are obtained by replacing p by a and g% by g\ everywhere. The flow generated by 
these equations can easily be seen to lead to a phase transition by plotting it in the gs — K p plane, as shown in Fig. 
|IV C| The flow for negative gs is analoguous, as can be inferred from Eqs. [72] Whilst the cosine term appears to 
be irrelevant for K p < 1 on dimensional grounds, it may actually be marginally relevant because K p is renormalized 
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FIG. 4: RG flow in the charge sector near K p — 1. For a massive spin sector, the diagonal in the right half is the separation 
line between the phases with leading sSC correlations (#3 flows to zero) and CDW order in the ground state (#3 is relevant). 



during the flow. This leads to the separatrix given by the straight lines in Fig. |IV C[ separating a phase where g% 



flows to strong coupling from a phase where it flows to zero. The system flows to strong coupling if 

|<73 1 > MK P - 1). (73) 

In the bosonization picture, the flow to strong coupling can be interpreted physically as the development of a gap in 
the charge boson spectrum. This follows from the fact that as g% becomes large, the cosine potential becomes steeper 
and steeper, so that the field <p p will get trapped in one of the minima, because the kinetic energy is no longer high 
enough to overcome the energy barriers separating the minima. Of course, one could object that tunneling between 
the minima via instantons might spoil this picture. However, it can be shown that finite energy instantons for the ID 
sine-Gordon model do not exist 25 . For large #3, the cosine may be expanded around the minimum, leading to a mass 
term (i.e. a gap) for </> p , 

Vf ^ 3 cos ^V^7n/>p^ « const. + ^ Vf ^ 3 ^ = const. + m(j) 2 p1 (74) 

where the mass, m, is defined by the above relation. 

Having found the low-energy behavior in the bosonic picture, the correlations in the fermionic picture can be 
discussed. Generically, one or more order parameters develop power law correlations. Due to the Mermin- Wagner 
theorem in one dimension continuous symmetries can not be broken (even at T = 0), so that no long range order 
develops if no discrete symmetry breaking is involved. Instead, one finds a scale-free state, a Luttinger liquid, that 
may be viewed as a system on the verge of various symmetry breaking instabilities. Typically, the degree of divergence 
of the static susceptibilities is not the same for the different instabilities, so that one instability is dominant. When 
coupled to a higher dimensional environment, the system would order according to its dominant instabilit^H. The 
susceptibilities can be expressed in terms of correlations of the corresponding order parameters, which in turn can be 
given in terms of the bosonic fields. Thus we will start with the bosonized expressions of the two order parameters 
that will be of interest later: s- wave superconductivity (sSC) and charge density wave (CDW). 

From the bosonization rules introduced above, it is easy to verify that the order parameters for CDW and sSC are 
given by 



OCDWW = e —^^* P cos (v^r) , (75) 
O sSC (x) = ^e- iv ^ cos (v^&t) • (76) 
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The correlation functions for these order parameters can be obtained from the correlators of the various vertex 
operators appearing in their definition. They read 

2ik F x , a .K p +K a 

<<W^CDW(0)> = ^ (") - 

(Ol sc (r)O sSC (0)) = ^{^) K ' 1+K \ (77) 

The behavior of the susceptibilities can be directly calculated from the order parameter correlations, leading to the 
singular behavior 

XCDw( k = 0,^) oc uF*+K--* (78) 

X sSC (fc = )W )a^" 1+ ^- 2 . (79) 

It should be noted that above formulae hold if both c/) a and (j) p are massless. Of course, K p and K a are renormalized, 
so actually one should first perform the RG calculation and afterwards use the renormalized values to calculate the 
singular behaviour of the susceptibilities. 

The field <j) a will be massive if \g\\ > 2i\(K a — 1) which is certainly the case for K a < 1, i.e. gu — g<i > 0. Later, 
in the fRG treatment, we will focus on this regime. Then the behavior of the spin-sector is modified as compared 
to the massless case. A massive <\> G implies that operators like cos^ acquire nonzero expectation values, or that 
their correlators are strongly suppressed, depending on their value at the potential minimum. For the attractive case, 
gi < 0, it turns out that all order parameter correlations apart from sSC and CDW are suppressed if </> a orders 24 . 
The correlations of Oqj^^jt and O s qq in this case are modified according to 

( CDWW CDW(0))«(^) KP 

<<?1SC WSC(0)> « © V , (80) 



and the same modification occurs for the static susceptibilities, 

XCDW x uKp ~ 2 

X sSC « " K;1 ~ 2 - (81) 

This obviously originates in the fact that a massive field can no longer fluctuate at arbitrarily low energies, so that 
(j) a can no longer influence the asymptotic behavior of the system. 

sSC will be the dominant instability as long as <p p remains massless, because in this regime K p > 1 always holds 
(where equality can only occur for g% = 0), so that X^SC a ^ wa y s diverges stronger than XCDW* 

What happens if (j) p becomes massive, too? The answer can be read off the bosonized form of the order parameters. 
O s qq depends on the dual field, 6 p . If <\> p orders, its dual field has to disorder, so that its correlations are exponentially 
suppressed^, and the sSC susceptibility is no longer singular. OQj^yy, on the other hand, acquires a finite value as it 
4>p. At half-filling, CDW only breaks Z2 which is a discrete group, so there are no low- lying modes that could destroy 
long-range order, so that at T = 0, CDW ordering occurs. 

To summarize, bosonization for our half-filled attractive model with massive spin sector predicts that the transition 



from massless to massive (f) p , tuned by crossing the separation line in Fig. IV C either by increasing \gs\ or decreasing 
K pi is accompanied by a transition from a phase with dominant sSC instability and subdominant CDW power-law 
correlations to a phase with CDW ordering and exponentially decaying sSC correlations. 



V. FUNCTIONAL RENORMALIZATION GROUP APPROACH 



In this section, the fRG with counterterms will be applied to the attractive Tomonaga-Luttinger model with 
interactions given by Eq. [3] at half-filling. We focus on the attractive case with massive spin sector, as this gives 
us the possibility to study the competition between charge-density- wave and superconducting ordering tendencies. 
The goal is to test the applicability of the fRG with symmetry-breaking counterterms to systems where there are 
competing correlations, and where the truncated flow equations are not exact. The bosonization phase diagram 
established in the last section will be the guideline for checking the quality of the results obtained. A word of caution 
is appropriate in the beginning: As shown above, the different phases of ID systems do not correspond to different 
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ordered states in general, but to critical states with one or several diverging static susceptibilities. In the currently 
used setup, the fRG cannot reproduce this critical behaviour. Instead, because order parameter fluctuations are not 
included, one arrives at the ordered state corresponding to the strongest instability. For this reason, numerical values 
like the order parameter should be understood as energy scales, and as indicator for the dominant correlation. The 
parameters found by fRG for the transition point from the superconductivity dominated phase to the charge density 
wave phase are used as a check for the correctness of the results. We see the main strength of the method tested here 
in two-dimensional applications, where long-range order in the ground state is more frequent than in one dimension. 
Hence the mischaracterization of power-law odered states as long-range orderd states is for us a tolerable defect which 
is clearly rooted in the nature of the approximation. 



A. Counterterms and Propagator 

As the model exhibits two different instabilities, we will need two counterterms to deal with them and to render 



the flows finite in the whole parameter range. These are the anomalous sSC self-energy already discussed in Sec. II B 



and the counterterm corresponding to a commensurate charge density wave (CDW). The order parameter for CDW 
symmetry breaking at half-filling, A^jyyy couples to the fermions via 



H CDW = A CDW E E C l,* c k+Q,< 



where Q = 2kp = 7r for half filling . In the language of chiral fermions, the relations 

Rk±Q = L k , L k ±Q — R k 

hold, so that Hqj^^ can be written as 

_ a J z?t r , _l r.t 



A CDW E E { R iM,° + L i* R k,*} • 



Z CDW 

Using similar relations, one finds for the pairing part H s $q that 



# S SC = A S SCE{4, T < 



Rt rt 



h.c.}. 



(82) 



(83) 



(84) 



(85) 



Next, we introduce a Nambu spinor notation for the fields, so that the free inverse propagator can be summarized 
in a 4 x 4 matrix structure. To this end, one must take care to arrange the fields in the spinors in a way such that 
^ and \P are independent from each other, because otherwise the mathematical structure of the action is changed 26 . 
The choice of spinors in this work is 




L_, L + , i?_ 



(86) 



where + is shorthand for (fc, j) and — means (— fc, j). The inverse free propagator including counterterms reads with 
these conventions 



Q 



i 




Ksc 







\ 






Q-k 





—Acdw 











Qk 


A sSC 




V 





—Acdw 


AsSC 


Q-k 


) 



(87) 



where 



Q_ k =iu - v F k, Q k = Q_ k = iuj + v F k 



represent the symmetric inverse free propagator, and Q is the transpose of Q. 

Introducing the flow parameter x £ [0? 1] f° r the interaction flow (as in Subsection III A), the full propagator is 

given by the standard formula Q = {^{Q — X^)) • Here E is again reserved for the flowing selfenergy. Herein we 
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take the counterterms Acdw and A s sc as initial conditions at x = 0- This exactly cancels the As added to Q for 
X = 1, so that in the end the model is unchanged. This leads to the expression 



X 



x 



(j 2 + (A sS c ~ X^ssc) 2 + (Acdw - X^cdw) 2 
I -icu + -e k A sS c - X^ssc Acdw ~ X^cdw 

^ssc ~ X^ssc -iu + e k 

A C dw - X^cdw —iw + e fc 

\ -Acdw + X^cdw ^ s sc - X^ssc 



(89) 





-Acdw + X^cdw 
A S( sc - X^sSC 
-iu - e k 



for the full propagator that occurs in the flow equations. One should note that the elements of the counter diagonal 
are all identically zero, implying that, apart from the symmetry breaking propagator entries for sSC and CDW, no 
new anomalous entries are created. The new propagator stemming from the CDW symmetry breaking correlates left- 
and right-movers of the same spin and will thus subsequently be called umklapp propagator. 



The single-scale propagator is obtained from the definitions of Q and Q above from its definition (cf. eq.(24), 



S = -QQG, 

where the product here is of course given by matrix multiplication. The entries with indices ij for which 

Qij = A x 

holds will be called S x with X = sSC, CDW. 



Of course, as already mentioned in section II B , the off-diagonal entries of the propagator may generate new vertices 
during the flow. In the case at hand, however, no new vertices are generated apart from the anomalous superconducting 



vertex introduced in section II B This vertex will be named g$ from now on, so that the full collection of interaction 
vertices is g u g 2 , #3, 9\\ , and g 5 . 



B. Flow Equations 



Using the definitions of the propagators and vertices of the above sections, one can directly derive the flow equations 
for the model. The result is 

£ a sc = \s a sc(9i+92-2g 6 ) (90) 

Scdw = ^Scdw(9i + 93 ~ 9\\) (91) 

for the order parameter flows (off-diagonal self-energies). Here ^sSC^CDw) denote 1°°P integrals over the corre- 
sponding matrix elements of the single scale propagator defined in the last section. 
For the interaction vertices the flow equations are given by 

gi = -( A + D )gl _ \ C g\ - Dg m - X -Cg\ - (C + D)g\ - gi (Cg 5 + Bg\\ + Ag 2 ) 

-92(Cg 5 + ^Dg 2 ) 

g 2 = -(A + D)g\ -\{B + C)g 2 3 - Cgj - C 9l g b - l -Ag\ - g 2 (Cg 5 + Dg 3 - 

Dg\\+ D 9l ) ~\{A + B + 2C)g 2 2 
9\\ = -\Bgl + (2D - X -A- X -B)g\ - Dg 3 g\ - X -Bg x - (Cg 5 + Dg 3 - Dg\\ + D 9l )g 2 

- l -(A + B + 2C)g 2 2 
g 3 = -Dgl - g 3 (B g]l + C gi + (B + C)g 2 ) - Dg ]igi - -Dgl 

£5 = ^{gl + g^-Cg^-Cgt-iA + D^+g^. (92) 
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The capital letters are abbreviations for the loops occuring on the RHS of the flow equations. They are explicitly 
given by 

A = ^d x T^2 u J dkQn(uj, /c)^33(o;, k) particle-hole loop 

B — \dxT^2uj J dkQn {jjj , k)Qss(— uj, —k) particle-particle loop 
C= ld x TJ2 (JJ Jdk(g 2 i(u; : k)) 2 sSC loop 

D = \d x T Z u J dk (g sl (<j, k)f CDW loop. 

The Matsubara sums appearing in the loops may be evaluated analytically, whereas the momentum integrals have to 
be evaluated numerically in general. However, at T = and with a linearized dispersion, all integrals can be obtained 
analytically, although they are quite complicated. 

Finally, following Gersch et alP^, we also calculate the free energy of the system, which is essentially determined by 
the flow of the zero-point 1PI vertex and the free energy of the initial, non-interacting, sytem. The equation can be 
written aeP 

n = -^Tr(XG), (93) 

where the intial condition Q(x = 0) = tto is calculated by integrating out the non-interacting fermions at x = 0- Here 
of course the symmetry breaking terms in the quadratic part of the action have to be taken into account. 



C. Self-Consistent Approach: Results Depend on the Choice of Counterterm 



When the truncation of the flow equations is not exact (as is the case for generic models), the results of the fRG 
with counterterms will in general no longer be counterterm independent. In fact, this is what happens for our example 
of the Tomonaga-Luttinger model. The problem occurs even if one investigates regions of the phase diagram where 
only one instability is dominant. This is shown in Fig. V C| for the case of dominant sSC instabiltiy for temperatures 
the critical temperature T c where symmetry breaking sets in. In other words, the selfenergies at the end of the flow 
become function of the counterterms A. Hence we write A). One also finds that for low temperatures not all 
choices of counterterm lead to convergent flows. 
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FIG. 5: Counterterm dependence of the sSC gap Ass — ^ss deep in the sSC regime (solid line). Also shown is the largest 
coupling gmax encountered during the flow (dashed line). 



Thus a method is needed to choose one value of the counterterm for which the results are expected to be optimal 
in some sense. We will follow an analogy to perturbation theory, or more precisly a method for determining optimal 
counterterms based on a variant of renormalized perturbation theory. This was used by Neumayr and Metzner^ 
to circumvent the occurence of divergencies in the renormalized perturbation expansion of many-fermion systems at 
T = 0. The problem there is that terms like 

GoTi • • • QoTiQo 

occur in the expansion which lead to divergencies at the Fermi surface because of the singular behaviour of the Green's 
functions. The strategy followed by Neumayr and Metzner is to choose a different starting point for the perturbation 
expansion by introducting a counterterm 5TL with 



Tio Tto + STl 



(94) 
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and to choose SH such that the self-energy vanishes everywhere on the Fermi surface, i.e. 



E(*0k=o = O. 



(95) 



This equation is solved iteratively to the desired order in perturbation theory, thus enabling one to perform calculations 
that remain finite as the zeroes of E cancel the poles of Qo at the Fermi surface. 

The physical idea behind this method is to choose an optimal starting point for the perturbative expansion, that is 
one where the non-interacting part already describes the system as well as possible. It is very natural to assume that 
the quality of a starting point - i.e. a choice of Ho - is measured by the amount that Ho is renormalized. If it is not 
renormalized at all, which is ensured for low energies by eq.(95), it is obviously already a good description of the full 
interacting model. 

The connection to the fRG approach to broken symmetries is established by noting that the method still works for 
symmetry breaking self-energies, such as superconductivity^. Furthermore, in the fRG scheme with counterterms 
one makes exactly the same division of the system into an interacting and a free part as in Eq. [94] As the truncation 
of the flow equations corresponds to a weak-coupling approximation, one is also interested in having a good starting 
point for the perturbation expansion, so that higher order corrections can be expected to be small. 



Based on these similarities, we propose to choose the counterterm A using the same condition as in eq.(95), which 
in fRG terms reads 



S( X = 1,A) = 0. 



(96) 



This equation is understood to hold for selfenergy appearing in the interacting Green's function as x^ with initial 
condition E(x = 0, A) = A. It does not correspond to a trivial situation, as the total or effective selfenergy, i.e. the 
difference to the non-interacting quadratic part, at x — 1 is given by E — A. 



Of course, Eq. [96] again can only be solved iteratively by inserting the fin al res ults for the negative total selfenergy 
-(E — A) at x = 1 as counterterm into a new run. However, as shown in Fig. V C , the iteration converges very quickly 



to a unique value. Thus it is possible to make a well-defined choice of counterterm out of a continuum of possibilities 
at moderate computational costs. In the case of the Tomonaga-Luttinger model considered here, typically less than 
ten iterations were needed for the method to converge. 




20 30 
iteration 



FIG. 6: Convergence of iterative solution of the self-consistency condition for g\ — —0.3, gz — —0.15, in units where vf = 1, a) 
in the sSC regime (#2 = 0), b) in the CDW regime (#2 = —0.5), and c) close to the transition (#2 = —0.25). 
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D. fRG Results 



As discussed before, for g 3 < 0, the bosonization description yields a quantum phase transition (QPT) in the 
K p — #3-plane between the power-law superconductor and the charge-density wave ground state. The question is 
wether the fRG scheme can reproduce this QPT. Hence we run the fRG as described above for a range of system 
parameters near the QPT. We use two criteria for determining the position of the phase transition: First, we use 
either the CDW or the sSC counterterm separately, and also calculate the free energy during the flow. We then 
choose the solution with lower energy (as done in Gersch et alP^) as the best approximation to the true behaviour. 



The resulting free energies are shown in Fig. VD for gs = —0.5 (in units where the lattice constant is unity and 
vp = 1)- We can read off a critical value of K p w 0.017 where the energetically favorable order parameter changes 
from pairing to CDW. Second, we can use both counterterms simultaneously, allowing thus for a coexistence of both 
orders, and say that for Acdw(x — 1) > ^sc(x = 1) the system is in the CDW phase and vice versa. In this case, 
the free energy for flows with order parameters allowed does only differ a lit tle from the curves for only one order 
parameter in a small coexistence region between 0.015 and 0.019 (see Fig. VD). The superconducting order parameter 
suppresses the CDW order on the large- if p -side, while the situation is reversed at small K p . This behavior of the 



order parameters versus K p is shown in Fig. VD Again we get a transition point of K p w 0.017. Hence, we find 



that the results obtained by these two ways to determine the correct counterterms (via order parameters and via free 
energies) coincide in the selection of the dominant order parameter. 

Next, as a further check of the reliability of the result, we use the fact that in the bosonization method the spin and 
charge sectors decouple, so that the phase diagram is independent of the value of g\ (as long as the spin sector is in the 
massive phase), and only depends on the combination gn — g2 entering K p . We find that even though this separation 
into different sectors is not manifest in the fermionic formulation, the position of the phase transition changes only 
slightly if g\ is changed with K p and gs fixed. This is shown in Fig. V D All this evidence taken together demonstrates 
that the fermionic fRG with counterterms is able to handle systems with more than one instability in a satisfactory 
way. 

On the back side it should be mentioned that near the transition between SC and CDW regime the fRG predicts a 
regime with coexisting order, where both order parameters are nonzero at the end of the flow. This regime is absent 
in the bosonization, where the incompressible CDW has exponentially decaying charge correlations. The fRG results 
do never give a gapless fermionic spectrum and replace the quantum critical point by a small coexistence region. 
While this difference is presumably due to the approximate treatment of the competition between the two ordering 
tendencies via static order only, it is certainly not so severe as to render the fRG approach useless. The ultimate nature 
of quantum phase transitions is often a subtle issue. Furthermore, as a more severe but also understandable difference 
of the fRG picture to the bosonization, the power-law superconductor is replaced by a long-range ordered symmetry- 
broken state. This was alluded already in the beginning of this section and is clearly due to our approximation in which 
only frequency-independent selfenergies are allowed, and where the simplified wave-vector- and frequency-structure 
of the vertices excludes a correct description of collective fluctuations. 

In Fig. 10 we show typical flows in the sSC and the CDW regime. In both regimes, all couplings remain moderate, 
such that the truncation error should remain acceptable. Couplings in the massless phase, where superconductivity 
dominates, become larger due to the presence of the Goldstone mode 13 . In an exact treatment one would actually 
observe a divergence of the couplings which reflects the fact that the Goldstone mode is massless. In the approximate 
treatment, this seems to be violated to some extent, i.e. the final sSC order parameter may be too large compared 
to the unknown exact result. In other words, the bias or explicit symmetry breaking of the initial conditions has not 
been removed entirely in the approximate flow. Of course one could try to vary the counterterm so as to recover 
the Goldstone mode. However it is unclear if this scheme would actually work, and moreover this alternative scheme 
for selecting the right counterterm could not be generalized to discrete symmetry breakings. We believe that for 
our purposes the imprecise treatment of the Goldstone mode is a quantitative problem that can be reduced by more 
sophisticated treatments of the frequency and momentum dependence of the vertex functions. Due to the qualitative 
correctness of the transition line between dominant CDW and superconducting correlations, the counterterm method 
in its present form should still be to determine the main picture. 



VI. CONCLUSIONS AND OUTLOOK 



Summarizing the above results, we have investigated the applicability of the interaction (flat cutoff) fRG flow scheme 
using counterterms to situations with symmetry breaking and competing orders. We have found that this method can 
handle problems where there is only on channel causing instabilities (the reduced BCS model) perfectly well. We then 
went on to apply the scheme to a system with competing instabilities, which led to several new insights. First of all, 
the results at the end of the flow are no longer independent of the choice of counterterm, so that one has to find a way 
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FIG. 7: Free energy for flows using only one counterterm (CDW or sSC) at a time (gray lines) and for flows with both 
counterterms (red line) in the vicinity of the phase transition. g% — —0.5 is kept fixed whilst K p is varied. 




FIG. 8: Behavior of the CDW and sSC order parameters in the vicinity of the phase transition, — —0.5 is kept fixed whilst 
K p is varied, for flows with both countertems included. 



to choose a preferred value of counterterm. We proposed to use solve a self-consistency equation iteratively, which 
turned out to give unique and sensible results in most cases. With this device at hand we investigated the quantum 
phase transition from a CDW ordered state to state with (power-law) superconducting order in a one-dimensional 
model with attractive fermions at half band filling. Two different ways to run the flows, either with only one type 
of symmetry breaking allowed, or with both types of symmetry breakings included gave the same location of the 
transition line from dominant CDW to dominant superconducting correlations, in agreement with the bosonization 
result. For most parameters, the coupling constants remain in the weak coupling range. 

Expected differences with respect to bosonization occurred in the long-range-ordered nature of the superconducting 
ground state in the fRG, and in the small coexistence region near the transition, smearing out the quantum criticality. 
Improvements in this regard would require a more precise treatment of collective fluctuations, which is beyond the 
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FIG. 9: Phase diagram in the K p — g% — plane for different values of g±. The bosonization result (solid line) is independent of 
<7i, while in the fermionic RG calculation there is a slight dependence on its value (colored points, g\ — 0.5, 1, 1.5, 2). 





FIG. 10: Typical fRG flows a) in the massless (SS) phase and b) in the massive (CDW) phase. The scale for the order 
parameters is on the left side, the scale for the couplings (multiplied by —1) on the right side. 



scope of this work. Nevertheless, the present method is able to predict correctly the type of leading ground state 
correlation and to produce a gapped renormalized single-particle spectrum. All fermionic modes can be integrated 
out without a severe run- away flow of the couplings. As a next step, the scheme should be applied to two-dimensional 
models. There it should allow one to compute, e.g., the order parameter magnitude and the gap structure around the 
Fermi surface of unconventional superconductors. 

We thank Manfred Salmhofer, Walter Metzner, Jutta Ortloff, T. Maurice Rice and Manfred Sigrist for useful 
discussions. This work was supported by the DFG research unit FOR723. 
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